Magnetic resonance from the interplay of frustration and superconductivity 
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Motivated by the iron-based superconductors, we develop a self-consistent electronic theory for 
the itinerant spin excitations in the regime of coexistence of the antiferromagnetic stripe order 
with wavevector Q 1 — (tv,0) and s + ~ superconductivity. The onset of superconductivity leads to 
the appearance of a magnetic resonance near the wavevector Q2 = (0, n) where magnetic order is 
absent. This resonance is isotropic in spin space, unlike the excitations near Qi where the magnetic 
Goldstone mode resides. We discuss several features which can be observed experimentally. 
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The idea of emergent symmetry in strongly corre- 
lated electron systems plays a prominent role in var- 
ious novel superconductors such as high-T c cuprates, 
heavy fermions, and organic compounds. The recent 
discovery of superconductivity in vicinity of the antifer- 
romagnetism in iron-based systems [l[ has opened up a 
new perspective in studying electronic instabilities and 
their symmetry properties. For example, it was argued 
initially that the close proximity of the antiferromag- 
netism (AF) and superconductivity (SC) in the iron- 
based superconductors could be understood within an 
emergent SO (6) symmetry of the underlying electronic 
Hamiltonian[2j, a result supported also by a renormaliza- 
tion group study. 3] This symmetry is closely connected 
to the so-called s H pairing state of the SC order param- 
eter. In addition, the frustrated AF state with ordering 
wavevector Q x = (n,0) or Q 2 = (0,7r) is related to an 
SU(2) x Z 2 symmetry where the Ising degree of freedom 
(Z 2 ) refers to the (tt,0) and (0, 7r) ordering vectors. This 
is analogous to behavior found in localized magnetic sys- 
tems, most simply in the frustrated J\ — J 2 square lattice 
model. 3 

As such enhanced symmetries are generically not pro- 
tected, fine details of the electronic structure lift this de- 
generacy, but only on a scale well below the Fermi energy 
(about 0.1 eV), and in a way which is "non-universal". 
Even only the selection of the AF stripe order out of the 
degenerate manifold of the magnetic structures associ- 
ated with two magnetic order parameters at wave vec- 
tors Qi and Q 2 and connected by the symmetry opera- 
tions of the 0(6) group represents an interesting problem. 
There, finite ellipticity and interactions between electron 
pockets remove the degeneracy and favor the observed 
stripe AF state. [j| Nevertheless, the competing states in 
the manifold remain nearby in energy^ @], resulting in 
collective excitations at low energy as a remnant of the 
enhanced symmetry. 

AF long range order naturally gives rise to a Goldstone 
mode. SC of the s H -type leads to a magnetic resonance 
mode instead. |8|, |9( It is an open question how these two 



rather distinct magnetic collective modes interfere in the 
case of coexisting AF + SC orders. Of course, the exis- 
tence of a gapless Goldstone mode (assuming continuous 
spin symmetry and singlet SC) is guaranteed based on 
general symmetry principles. It seems plausible to spec- 
ulate that that the resonance mode is simply absorbed by 
the Goldstone mode. Here we show that the enhanced 
Z 2 -symmetry and the existence of two potential order- 
ing vectors does allow for a strong resonance mode at 
Q2 = (0,7r) in the coexistence phase if the magnetic or- 
der corresponds to Qi = (tt, 0) and vice versa. This ef- 
fect is peculiar to itinerant magnetic systems with nearly 
degenerate ground states and demonstrates that the col- 
lective behavior of various excitations may uncover the 
hidden symmetry of the electronic system. 

The electronic structure of the iron-based supercon- 
ductors can be approximated by the circular hole pock- 
ets centered around the T = (0, 0) point of the Brillouin 
Zone (BZ) and elliptic electron pockets centered around 
the (tt,0) and (0,7r) point of the BZ (with 1 Fe ion per 
unit cell). For our purposes, it suffices to consider the 
simple two-band model[lO| 

H = E {^)cLc ka + e/(k) jt/ kCT } + 

k,o- 

E {wi4 +qCT /^,_ qCT ,jWc kCT + 

y (/k+ qCT /k'- qff ' c k' CT 'Ck CT + /i.e.) } (i) 

where e c (k) = 2.0i (cosfc^ + cosk y ) +e c + fi and e/(k) = 
2.0t (cosk x ■ cosfcj,) — t' (cosk x + cosk y ) + e; + fj, refer 
to the hole (c) and electron (/) bands. We further set 
t = 1.0, H = 0.7t, e f = 1.3t, e c = -3.3t, and \x = 0.03t 
with the resulting Fermi surface shown in FigfTfa). The 
values of the hopping integrals slightly vary from those of 
Ref. [l(| to include finite doping and ellipticity. u\ and U3 
refer to the Fermi-liquid like interactions whichgive rise 
to the AF and SC order with symmetry. [3j In the 
following we concentrate on the coexistence phase 11 - IH 
in which s" 1 superconductivity coexists with antifcrro- 
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magnetic order in a rather broad range of parameters. 

Method. We decouple the interaction part of the 
Hamiltonian with respect to the AF and SC order, as- 
suming T c <TV- We first apply the unitary trans- 
formations with respect to the AF state and introduce 
the new quasiparticle operators Qk.a, Ac, a- for the two 
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cL/ka ) sqn(cr) with u S dw =ui+u s 



is the AF order parameter with Qi ordering momentum 
and = ec ( k )^ e ^ k ) , i n the AF state with ordering at 
Qi = (?r, 0) the system remains a metal as the pocket at 
(0, 7r) remains intact. In addition, because of the finite 
doping and cllipticity of the electron pockets, tiny pock- 
ets of electron and hole character remain around the T 
point of the BZ at intermediate values of the AF order 
as shown in FigfTJb). 




FIG. 1. (color online) Calculated Fermi surface topology in 
the normal (a) state and AF (b) state with Qi = (n, 0) order- 
ing momentum. The + and — signs refer to the phase of the 
superconducting s +_ -wave order parameter in the coexistence 
phase. The arrows refer to the AF ordering wave vector and 
the wavevector of the incommensurate resonance excitations 
below T c . 



In the next step we rewrite the SC pairing interaction 
in the s" 1 channel with the unitary AF transformation, 
and subsequently perform a mean-field (MF) decoupling 
in the particle-particle channel, keeping only anomalous 
expectation values of (oJoJ) and (/J'/?'). The resulting 
MF Hamiltonian is then diagonalized by two independent 
Bogolyubov transformations for the a and /3 band, yield- 



ing the new energy dispersions £1^ 



(A^ 2 



where 7 = a, (3. The SC gaps in the s H channel, A^ ar< 
determined self-consistently from two coupled gap equa- 
tions 



A, 



where Ai = 
A 2 = 



ulA 2 
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The folding of the electron 



2k B T 

and hole bands with opposite sign of the supercon- 
ducting gap [lH at Qi does not produce nodal lines 
at the intersection points, [lj] Finally, Eqs.([2]) have 
to be supplemented by the self-consistent equation 
for the AF order itself in the coexistence phase. 1 = 



7W 



mi) - \] } 



The use of the sequential transformations is valid when 
the magnitude of the SC gaps is small compared to the 
SDW gap. However, more generally it is possible to 
show that the enhanced symmetry of the problem makes 
this approach valid for arbitrary ratios |A,-|/W as long 
as the SC gaps of the two bands are similar. 

To obtain the spin susceptibility we employ a general- 
ized random phase approximation approach in the multi- 
band casep, [l(J UM- in terms of Green's functions, the 
dynamical susceptibility tensor for the longitudinal, zz, 
and the transverse, H — , components is defined as 



St 



r(q^) = -^££[G: 



^(»W n )G^^(iu; n + ifi)T 



p.CT 

F;%(iu n )Fl a + ^ l( ,,(iun + in)] , 

where m = zz, H — and a = a' (a = a') 
for the longitudinal (transverse) component. Here, 

G^(iw n ) = -fgdriTrS^iT^me™^, and 

F pn(iUn) = - Jo dT(T T s pt (T)t pl (0))e tw " T . Then, the 
spin susceptibility is obtained via a Dyson equation 

\Xba\RPA = Xba + iXb'a'} RPA ^ c'% Xba • 

In the following we evaluate numerically the gap equa- 
tions for the AF state with Qi = (tt, 0) ordering wave vec- 
tor and s H SC self-consistently at T=10K. We choose 
the following interaction parameters that allow for a self- 
consistent solution of the AF and SC gaps at this temper- 
ature: usdw ~ 4.9t and M3 « 5.5t. Correspondingly we 
obtain a AF gap W = 0.097t and SC gaps Ai = 0.0m, 
A 2 = 0.016i. We note that self-consistency is crucial to 
obtain the correct spectrum of the spin excitations in the 
AF and AF+SC state. 

Results. To build up the picture for the coexistence 
phase, we consider the situation of the pure AF and s H 
SC phases. In particular, setting SC gaps Ai = A 2 = 
we find W = 0.136t, significantly larger than the AF gap 
coexisting with SC because both orders compete for the 
same FS and suppress each other. In FigfSJa) we show 
the imaginary part of the transverse spin susceptibility 
around the magnetic ordering wave vector in the AF 
state only. The dispersion of the spin waves is strongly 
anisotropic in q x (rj < 0) and q y (r) > 0) direction due 
to ellipticity of the electron pockets involved, well visible 
up to an energy of 2V(^.0 Despite the mismatch between 
the sizes of electron and hole pockets separated by Qi, 
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FIG. 2. (color online) The (q, Q) intensity map (in states/t) of 
the imaginary part of the transverse spin susceptibility around 
the AF wave vector Qi = (n, 0) along a (ir, 0.157r) — !> (n, 0) — > 
(1.157r, 0) cut of the BZ, parametrized by -q, in the pure AF 
(a) and pure s" 1 superconducting (b) states. The inset in (a) 
shows the imaginary part of the transverse spin susceptibility 
around Q2 = (0, n). 



the spin waves are commensurate. 

In contrast to the gapless behavior of the transverse 
spin excitations around the ordering vector Qi, they are 
gapped near Q 2 as shown in the inset of FigHJ This gap, 
SI2 « 0. li, is a consequence of the ellipticity of the elec- 
tron pockets which lowers the symmetry of the model and 
favors the stripe magnetic structure over magnetic orders 
involving both momenta, Qi and Q2-|5, 6] For decreasing 
W, the spin excitations are still peaked around £7 2 but 
they are not necessarily gapless due to the onset of the 
particle-hole continuum, characterized by Q p h ~ W. The 
latter originates from the scattering between tiny rem- 
nant FS pockets and the electron Fermi surface centered 
around (Q,7r) as depicted in FigfTJb). For Vt p h << ^2, 
the excitations around Q2 are thus overdamped param- 
agnons. This is further supported by the fact that the 
value of Q2 is the same in the longitudinal, Imx zz , and 
the transverse, Imx" 1 , part of the spin susceptibility. 
This contrasts to the behavior near where the spin 
rotational symmetry is explicitly broken, i.e. the Gold- 
stone mode forms in Imx" 1 and the longitudinal excita- 
tions are gapped by twice the AF gap energy, 2W. 

Next we set the AF gap W = and solve the mean-field 
gaps for the superconducting s" 1 state self-consistently 
obtaining Ai « O.Oli and A 2 « 0.04t. We note that this 
limit has to be taken with care as our formalism initially 
relies on the assumption that W S> Aj or similar magni- 
tudes of the SC gaps on the electron and hole pockets. 

In addition, the correct behavior of the superconduct- 
ing gap for W — > has to be taken. The results are 
shown in Fig[2jb). Due to the fact that Ak = — Ak+Q AF 
in the s H superconductor, a dispersing resonance is seen 
at nonzero energy below twice the SC gap values. [3, 0] 
However, in contrast to the pure AF, the position of the 
resonance is at the incommensurate momentum, Q res ~ 



(1.057T, 0). The incommensuration is induced by the spe- 
cific shape of the FS with different sizes of the electron 
and hole pockets. [3] 

Finally, in Figj3] we present the results for coexisting 
AF+SC orders. In Fig|3ja) we show the evolution of the 
transverse spin excitations around Qi. Due to the finite 
magnetization, well-defined anisotropic spin waves exist 
around the ordering momentum similar to the pure AF 
state. This is expected as the Goldstone mode in the AF 
phase reflects the breaking of the spin-rotational sym- 
metry of the system and as such is not affected by the 
presence of the additional spin-singlet superconducting 
state. At the same time, and in contrast to the pure 
AF state, spin waves are damped away from Qi at en- 
ergies well below twice the AF gap 2W. In particular, 
additional damping starts at energies of twice the su- 
perconducting gap. This new feature of the coexistence 
phase arises due to the renormalization of the particle- 
hole (p-h) continuum. In the AF, the p-h response is 
determined by the transition between a— and /3— bands 
as well as the intraband transitions. The latter are sup- 
pressed by the vanishing AF matrix elements around the 
ordering wave vector. This also guarantees the stability 
of the Goldstone mode even in the situation when the 
AF gap does not gap the entire FS and remnant elec- 
tron and hole pockets are still present. The damping of 
the spin waves arises from interband transitions. They 
are gapped exactly by 2W at Qi; however, this gap be- 
comes progressively smaller away from it. In the coexis- 
tence phase, the interband excitations at Qi are gapped 
by the sum of twice SC and AF gaps. Away from Qi 
the effect of the AF quickly vanishes. However, the p- 
h continuum is still gapped by approximately twice the 
SC gap. Above this energy it experiences a discontinu- 
ous jump due to sgn(A£) = — sgn(A^). As a result, the 
spin wave excitations are damped at energies associated 
with twice the SC gap as seen from Figj3ja). Note that 
the 'stronger' damping of the spin waves in the coexis- 
tence phase [FiglSJa)] originates from the smaller value 
of the AF gap in the coexistence phase. The discontinu- 
ous jump of the particle- hole continuum at approximately 
2Ao is absent in the pure AF state. 

In Fig[3jb) we show the spin excitations in the coexis- 
tence phase around Q2. As discussed above, in the pure 
AF phase the excitations around Q2 are paramagnons 
centered at frequency f?2 and almost isotropic in spin 
space. In the coexistence phase, however, the SC gap 
modifies the behavior of the p-h continuum around Q2 
which was gapless in the AF only case when tiny hole 
pockets centered around (0, 0) were still present. The SC 
gaps the p-h continuum up to 2Ao- In addition, due to 
the change of sign of the gaps between the pockets, the 
p-h continuum experiences a discontinuous jump at this 
frequency and the necessary condition for the formation 
of the excitonic resonance at energies below 2Ao forms. It 
is remarkable that this resonance exist at energies which 
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FIG. 3. (color online) The (q, fi) intensity map (in states/t) 
of the imaginary part of the transverse spin susceptibility 
around the AFM wave vector Qi = (tt, 0) along (7r,0.157r) — > 
(tt, 0) — ¥ (1.157r,0) direction of the BZ (a) and around the 
Q 2 = (0,tt) along (0, 1.15tt) (0, tt) (0.157T, 7r) direction 
in the coexistence SC+AF phase. 



are much lower than the energy position of the original 
mode at 0,2 in the AF state. 

The incommensuration of the resonance in the coex- 
istence phase is closely related to the reconstructed FS 
for which hole pocket states that would be connected by 
Q 2 to an electron pocket are gapped for sufficiently large 
magnetic order. Note also that contrary to the ordering 
momentum, the excitations around Q2 remain isotropic 
in spin space, i.e. the excitations in the longitudinal and 
transverse channel show qualitatively similar behavior. 

Overall the results in the coexistence region show the 
peculiarity of the spin excitations in the ferropnictides: 
while around the ordering momentum Qi the excitations 
are determined by the presence of the AF order, the ex- 
citations around Q2 are dominated by the SC and the 
formation of the incommensurate spin resonance. The 
formation of the resonance effectively lowers the energy 
splitting between the different magnetic states in the co- 
existence region. Correspondingly the system tends to 
increase its magnetic degeneracy in the presence of s" 1 
superconductivity. As a consequence, the collective exci- 
tations associated with the presence of the incommensu- 
rate spin resonance close to Q2 may occur in the coex- 
istence region, a phenomenon which would be absent in 
the pure AF state. This possibility would be interesting 
to check experimentally. 

A direct probe to see the excitations around Qi and 
Q2 separately in INS experiments is to look at untwinned 
crystals. However, existing experiments on twinned crys- 
tals already indicate that both spin wave and spin reso- 
nance are present in the coexistence phase. Ref.[17] found 
that spin excitations in the coexistence AF+SC phase 
show strong anisotropy along the c-crystallographic di- 
rection, absent in the pure SC state. The difference in 
energy of the spin excitations at (rr,TT,q z ) is 3 — 4meV 
between q z = and q z — tt. At the same time, it 
was argued [18[ that due to the peculiar crystallographic 
structure of 122 compounds, the excitations at (tt, tt, 0) 



and (tt, tt, tt) in the folded BZ zone correspond to exci- 
tations at (7r, 0,0) and (0,7r, tt) wave vectors in the un- 
folded BZ, respectively. In other words, while excita- 
tions around (tt, tt, 0) refer to the true spin waves at Qi 
(gapped in the experimental system by spin-orbit cou- 
pling we have omitted) the excitations around (tt, tt, tt) 
originate from the incommensurate resonance at Q2 . The 
crucial test to see whether this is true is to study the spin 
rotational symmetry of these excitations and their possi- 
ble incommensurability. 

We have analyzed spin excitations in the coexistence 
phase of iron-based superconductors in which both AF 
and s" 1 SC coexist as a model of an approximate emer- 
gent symmetry in an itinerant magnet. We find that the 
excitations around the 'failed' ordering wavevector Q2 
become resonant at energies VL res < 2Aq. In addition, 
the longitudinal component of the spin susceptibility dif- 
fers strongly near Qi and Q2. Around Qi, where the 
transverse excitations are gapless Goldstone modes, lon- 
gitudinal excitations are gapped by 2W on account of the 
presence of AF order and are only weakly affected by the 
SC; while at Q2, excitations are isotropic with a gap now 
set by the superconductivity, 2A . 
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